Exploring behaviors of partonic matter forming in early stage of pp and 
nucleus-nucleus collisions at ultra-relativistic energies 
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The parton yield, (pseudo)rapidity distribution, and transverse momentum distribution in par- 
tonic matter assumed forming in the early stage of pp and nucleus-nucleus collisions at RHIC energy 
(v / sjv"iv=200 GeV) and LHC energy (y / sjvjv=5.5 TeV for nucleus-nucleus, 5.5 and/or 14 TeV for pp) 
are comparatively investigated with parton and hadron cascade model PACIAE. It turned out that 
the different parton and anti-parton spectra approach to be similar with reaction energy increasing 
from RHIC to LHC. We have argued that if the partonic matter forming in Au+Au collisions at 
RHIC energy is strongly interacting quark-gluon plasma, the one forming in Pb+Pb collisions at 
LHC energy might approach the real (free) quark-gluon plasma. 

PACS numbers: 25.75.-q, 24.85. +p, 24.10.Lx 
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The coalescence (recombination) models [T|, 0, IE 0, H, 
Q have achieved great successes in relativistic heavy-ion 
collisions. In these models, different to the string frag- 
mentation model, the relevant degrees of freedom are not 
"free" quarks but massive (dressed) quarks. The gluons 
are assumed to convert to quark pairs. It is always as- 
sumed in these models that the spectrum of coalesced 
hadron is proportional to the spectrum product of the 
coalescing partons and/or anti-partons. Here the parton 
refers to d, u, s, ... quarks, and gluons, anti-parton refers 
to d, u, s, ... . The coalescence (recombination) models 
require the knowledge of parton pt and/or 77 distribu- 
tions as inputs 0. Similarly, the perturbative quantum 
chromo-dynamics (pQCD) calculations [j| Q for hadron 
production in the relativistic elementary and/or nuclear 
collisions need the intrinsic transverse momentum distri- 
bution of parton in the nucleon as well. 

Based on the further assumption of different quark and 
antiquark have same spectrum, coalescence (recombina- 
tion) models predicted that the elliptic flow parameter 
v-2 follows a quark-number (n q ) scaling in i^Cpt/?^) vs. 
PT/n q for most final state hadrons in the intermediate 
Pt region H, 0, E3, EH • Later, this quark- number scaling 
prediction has been proved experimentally [H, UH . This 
V2 quark-number scaling observation relates the spec- 
trum of final state hadron to the spectrum of initial state 
parton directly. 

Recently, the concept of effective constituent quark 
has been introduced to connect with the quark coa- 
lescence (recombination) picture [Tj, They as- 
sume that the hadron's pt is composed of its effec- 
tive constituent quark's pr/n (n is the number of ef- 
fective constituent quarks in hadron) and the different 
quark and anti-quark have same px distribution. Then 
they extracted the pt distribution of effective u (d) 
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and s quarks from the ratio of ^(pt/3)/0(pt/2) and 
£1{pt I '3) / '4>{pt 1 '2) , respectively, in Au+Au collisions at 
7i^7=200 GeV. 

As mentioned in [lj| recently that, "... there has been 
an effort to fully reconstructed jets in heavy-ion colli- 
sions in order to provide a direct measurement of the 
partonic kinematics, independently of the fragmentation 
process (quenched or unquenched)." In fact, experimen- 
tally extracting the partonic observables from final state 
hadrons with reconstruction method have already be- 
come interesting [l6l [l7j . So it is important to investi- 
gate the properties (yield, px distribution, and 77 distri- 
bution) and explore their different behaviors in pp and 
nucleus- nucleus collisions at RHIC energy (y / J/\w=200 
GeV) and LHC energy ( v / J/vtv=5.5 TeV for nucleus- 
nucleus, 5.5 and/or 14 TeV for pp) theoretically. These 
studies might shed light on the way toward QGP. 

In this paper the parton and hadron cascade model, 
PACIAE is used to investigate systematically the 
yield, (pseudo)rapidity distribution, and transverse mo- 
mentum distribution of partons and anti-partons in the 
partonic matter assumed forming in the early stage of 
pp and nucleus- nucleus collisions at RHIC and LHC en- 
ergies. The discrimination, in the yields, rj distribu- 
tions, and pt distributions, between different parton and 
anti-parton, RHIC and LHC energies, as well as pp and 
nucleus-nucleus collisions is discussed. We have turned 
out that the different parton and anti-parton spectra 
approach to be similar with reaction energy increasing 
from RHIC to LHC. It is argued that if the partonic 
matter forming in the early stage of Au+Au collisions 
at RHIC energy is a strongly interacting QGP (sQGP) 
0, 113, HH HI the one forming in Pb+Pb collisions 
at LHC energy might approach to the real (free) QGP 
(fQGP). 

PACIAE is a parton and hadron cascade model [l8j 
based on PYTHIA [23]. PYTHIA is a model for hadron- 
hadron collisions. The PACIAE model is composed of 
four stages of the parton initialization, parton evolu- 



tion (rescattering), hadronization, and hadron evolution 
(rescattering) . 

A nucleon-nucleon (NN) collision in the PYTHIA 
(PACIAE) model is decomposed into parton-parton col- 
lisions. The hard and soft parton-parton collisions are 
described, respectively, by the lowest-leading-order (LO) 
pQCD parton-parton cross section [24j and an empir- 
ical method. The semihard (between hard and soft) 
QCD 2^2 processes are involved as well. Because the 
initial- and final-state QCD radiation are considered, the 
PYTHIA (PACIAE) model generates a multijet event for 
a NN collision. This is followed by Lund and/or Inde- 
pendent Fragmentation model in the PYTHIA model, so 
one obtains a hadronic state for a NN (hadron-hadron, 
hh) collision. Since above fragmentation is switched-off 
in the PACIAE model, so one obtains a multijet event 
(composed of quarks, anti-quarks, and gluons) instead. 

A nucleus-nucleus collision in the PACIAE model is 
decomposed into NN collisions according to the colli- 
sion geometry. The nucleons in a colliding nucleus are 
arranged randomly in coordinate space according to the 
Woods-Saxon distribution (for radius r) and 47r isotropic 
distribution (orientation). We assume p x = p y = and 
p z equals the beam momentum for every colliding nu- 
cleon. Assuming straight line trajectory for nucleons we 
can calculate the collision time for each AW collision pair 
provided that the closest approaching distance between 
two colliding nucleons is less than or equal to y^f/f 
(a^ refers to the total cross section of NN collision). 
If every JVJV collision is performed with the method in 
previous paragraph until the collision pair is exhausted 
, we obtain a initial partonic state for a nucleus-nucleus 
collision. 



TABLE I: Composition of the partonic matter forming in 
pp and 0-5% most central nucleus-nucleus collisions at RHIC 
and LHC energies. 



PP 



Au+Au Pb+Pb 



Energy (GeV) 
d 
d 
a 
u 
s 
s 
9 



200 14000 200 



5500 



2.09 
0.237 
3.98 
0.234 
0.123 
0.106 
4.82 



5.25 
3.47 
7.08 
3.47 
1.19 
1.17 
39.8 



330 
27.3 
294 
27.5 
13.5 
12.7 
516 



244 
112 
226 
112 
39.6 
39.5 
1351 



In the parton evolution (partonic rescatterin g) s tage 
the 2^2 LO-pQCD differential cross sections [24j are 
employed. For a subprocess ij — ► kl the differential cross 
section reads 



= A 

dt s 
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FIG. 1: Pseudorapidity distribution of parton and anti- 
parton in pp collisions at -^5=14000 and 200 GeV. The solid 
and open circles are for partons in pp collisions at LHC and 
RHIC energies, respectively. Similarly the solid and dashed 
curves are for anti-partons. Panel (a), (b), (c), and (d) is for 
d and d, u and u, s and s, and g, respectively. The dashed 
curve has scaled by 5 in panels (a) and (b). 



where the K factor is introduced for the higher order 
corrections and the nonperturbative correction as usual. 
Take the process q\q2 —* 9i<?2 as an example one has 



E 



It can be regularized as 



E 



9 
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q\qi^>q\qi 
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(2) 



(3) 



by introducing the parton colour screen mass, /^=0.63 
GeV. In above equation s, t, and u refer to the Man- 
delstam variables and a s = 0.47 stands for the running 
coupling constant. The total cross section of the parton 
collision i + j is then 



a tJ (s) 



E 

k,l 



dt 



da. 



ij — >kl 



dt 



(4) 



With the total and differential cross sections above the 
parton evolution (rescattering) can be simulated by the 
Monte Carlo method until the parton-parton collision is 
ceased (partonic freeze-out). 

The hadronization follows parton evolution. In the 
PACIAE model partons can be hadronized by the string 
fragmentation scheme or the coalescence picture. As all 
the simulations are ended after parton rescattering in 
this paper, we do not describe the hadronization in detail 
but refer to Ref. 

After hadronization we obtain a hadron list composed 
of the spatial and momentum coordinates of all hadrons 
for a nucleus-nucleus collision. Similar to the parton 
rescattering above, one calculates the collision time for 
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each hh collision pair. Then each hh collision is per- 
formed with usual two-body collision method [25| until 
the collision pairs are exhausted (hadronic freeze-out). 

As we aim at the physics rather than reproducing the 
experimental data, so in the calculations default values 
given in the PYTHIA model are adopted for model pa- 
rameters except K—3 is assumed. The simulations are 
all ended at parton evolution (rescattering) stage. 
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FIG. 2: Transverse momentum distribution of parton and 
anti-parton in pp collisions at ^=14000 and 200 GeV. The 
solid and open circles are for partons in pp collisions at LHC 
and RHIC energies, respectively. Similarly the solid and 
dashed curves are for anti-partons. Panel (a), (b), (c), and 
(d) is for d and d, u and u, s and s, and g, respectively. 

Four white papers [H, HH HH have explained that 
the partonic matter (quark gluon matter, QGM) formed 
in the early stage of Au+Au collisions at v / saT/v=200 
GeV is sQGP. Suppose the QGM formed in the early 
stage of Pb+Pb collisions at v /s]vjv=5.5 TeV is a step to- 
ward fQGP. Then one could make conjectures for fQGP: 

• According to the momentum fraction carried by 
the gluons in proton total momentum is 0.5 [261 ] 
and inspired by the gluon saturation concept in 
the colour glass condensate model [27| , we suppose 
the above fraction in Pb+Pb collision at LHC en- 
ergy is lager than that in Au+Au collision at RHIC 
energy, because of the more gluon condensate in 
former reaction than in latter one. This fraction 
factor becomes larger than 0.5 with the QGM ap- 
proaching fQGP. 

• In the Lund string fragmentation regime, the fla- 
vor selection of u : d : s : c... s» 1 : 1 : 0.3 : 
10 -11 ... and uu : dd : ss, ... » 1 : 1 : 0.3 ... 
[23| is introduced in hadron production based on 
the mass (dressed mass) effect. As this effect de- 
creases with momentum (reaction energy /temper- 
ature) increasing (28|, it is reasonable to assume 
that a flavor balance of u : d : s : c... w 1 : 1 : 7^ : 



7 C ... and uu : dd : ss, ... w 1 : 1 : 7 S ■■■ (7s > 0.3 
and 7 C > 10 -11 for instance) is expected in fQGP. 

• It is well known the anti-particle to particle ra- 
tio approaches to one with energy (temperature) 
increasing, the ratio of anti-parton to parton is as- 
sumed approaching to one in fQGP. 

The transport model results in this paper support above 
three conjectures. 

In Tab. |T] we give the parton chemical composition 
of partonic matter assumed forming in the early stage 
of pp and nucleus- nucleus collisions at RHIC and LHC 
energies. One sees in this table that: 

• Assuming the sum of yields in each column is the 
total yield of partons in the corresponding reaction. 
The fraction of gluons in total yield is about 0.420 
and 0.422 in pp and Au+Au collisions at RHIC en- 
ergy while those are about 0.648 and 0.636 in pp 
and Pb+Pb collisions at LHC energy. Similarly 
the momentum fraction carried by gluons in to- 
tal momentum is about 0.165 and 0.179 in pp and 
Au+Au collisions at RHIC energy while those are 
about 0.420 and 0.415 in pp and Pb+Pb collisions 
at LHC energy. These results indicate that the 
QGM forming in the early stage of Pb+Pb col- 
lisions at LHC energy is closer to fQGP than in 
Au+Au collisions at RHIC energy. 

• In both pp and nucleus- nucleus collisions the yield 
of d, u, s, s, and g increases with reaction energy 
increasing from RHIC to LHC much stronger than 
d and u. As the yield of d, u, s, s, and g in Pb+Pb 
collision at LHC energy is a few times larger than 
the corresponding one in Au+Au collision at RHIC 
energy, the yield of d and u in Pb+Pb collision is 
even somewhat less than that in Au+Au collision 
in order to meet with the above three conjectures. 

• d and u yields are about a magnitude larger than 
d and u, respectively, in pp and Au+Au collisions 
at RHIC energy. But that yield difference drops 
to two dramatically in pp and Pb+Pb collisions at 
LHC energy. This is consistent with the conjecture 
of the ratio of anti-parton to parton approaches one 
in fQGP. 



In Fig. [T] we give the pseudorapidity distributions of 
partons and anti-partons in pp collisions at y / s=14000 
and 200 GeV. One sees in panels (a) and (b) that the 
fragmentation peaks survive and midrapidity valley ap- 
pears in the 77 distribution of d and u quarks in pp colli- 
sions at RHIC energy. But the fragmentation peaks ap- 
proach to disappear and midrapidity valley fills up with 
energy increasing from RHIC to LHC. However the frag- 
mentation peaks and midrapidity valley are not shown in 
the r\ distribution of d and u quarks. The 77 distribution 



4 



- d: circles; d: curves 
soiio. ro+rO. 1 e V 
open: Au+Ali 200 GeV 


- u(u) 


°o °° 

- Ca) 

o 
. ° o°. ' • 
•<_! * *' *" * * * 


O q O (b) 

: 00 o ° 








1 | 1 1 1 | 1 1 1 | 1 H | 1 1 1 | 1 1 1 | 1 1 1 | 

s(5) 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


" (c) /^**\' 


(d) 

" cF°o ' 




• o o • 
• o • 


m D \ \ 





- 120 

- SO 
40 




-4 4 
1 



-4 4 
1 



FIG. 3: The same as Fig.l7]but for 0-5% most central Pb+Pb 
and Au+Au collisions at y / 5jvjv=5500 and 200 GeV instead 
of pp collisions. The dashed curve has scaled by 2.5 in panels 
(a) and (b). 



of u quark is higher (RHIC energy) or wider (LHC en- 
ergy) than d quark because proton is composed of uud. 
However the rj distribution of u quark is nearly the same 
as d in pp collisions at both RHIC and LHC energies. 
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FIG. 4: The same as Fig.[2]but for 0-5% most central Pb+Pb 
and Au+Au collisions at v / 5jvjv=5500 and 200 GeV instead 
of pp collisions. 

Globally speaking, we know from Fig. Q] that the r\ 
distributions of partons and anti-partons in pp collisions 
at LHC energy is more similar with each other than pp 
collisions at RHIC energy. So it is reasonable to suppose 
that if QGM is not possible forming in the early stage 
of pp collision at RHIC energy, it might be possible at 
LHC or higher energies. This idea is consistent with the 
knowledge of the nuclear system is able to deconfine by 
either highly heating or strongly compressing [29|. 

The transverse momentum distributions of partons 
and anti-partons in pp collisions at RHIC and LHC en- 
ergies are given in Fig. O We see in this figure that: 



• The transverse momentum distribution of d (it) 
quark is different from d (u), while the px distri- 
bution of s and s is nearly the same in pp collisions 
at RHIC energy. 

• In pp collisions at LHC energy the pr distribution 
of d (u) is the same as d (u) in the pr > 2 GeV/c 
region. But s and s have the same px distribution 
completely. 

• The slope parameter (apparent temperature) in 
the pt distribution of different parton and anti- 
parton in pp collisions at RHIC energy is larger 
(lower) than the corresponding one in pp collisions 
at LHC energy. 

These results indicate again that the different parton 
and anti-parton px distributions approach similar with 
reaction energy increasing from RHIC to LHC. So one 
can not rule out the possibility of QGM forming in pp 
collisions at LHC or higher energy. 

Similar to Fig. [1] we give the 77 distributions of par- 
tons and anti-partons in 0-5% most central Pb+Pb and 
Au+Au collisions at ^/sj^ =5500 and 200 GeV in Fig. [31 
respectively. A discussion for Fig.[3]parallel to Fig.Q]can 
be drawn. It has only to emphasize again that the shapes 
of partons and anti-partons rj distributions are more sim- 
ilar with each other in Pb+Pb collisions at LHC energy 
than in Au+Au collisions at RHIC energy in globally 
speaking. Therefore if the partonic matter forming in 
the early stage of Au+Au collisions at RHIC energy is 
sQGP the one forming in Pb+Pb collisions at LHC en- 
ergy might move forward to fQGP. 

Figure 2J similar to the Fig. [51 gives the pt distribu- 
tions of partons and anti-partons in 0-5% most central 
Pb+Pb and Au+Au collisions at LHC and RHIC ener- 
gies, respectively. In this figure one sees the same fea- 
tures as in Fig. [2j The pt distribution of d (u) quark is 
different from d (u) , while the px distribution of s and s 
is the same in Au+Au collisions at RHIC energy. How- 
ever in Pb+Pb collisions at LHC energy the different 
parton has the same px distribution as corresponding 
anti-parton (for d and u quarks in the pr > 2 GeV/c 
region only). The slope parameters (apparent tempera- 
tures) in the pt distributions of different parton and anti- 
parton in Au+Au collisions at RHIC energy are larger 
(lower) than the corresponding ones in Pb+Pb collisions 
at LHC energy. Thus if the partonic matter forming 
in the early stage of Au+Au collisions at RHIC energy 
is sQGP the one forming in Pb+Pb collisions at LHC 
energy might approach fQGP. 



In Fig. [5] we compare the different parton and anti- 
parton 77 distributions in pp collisions at ^/s=5500 
and 200 GeV with the corresponding ones in Pb+Pb 
and Au+Au collisions at the same cms energy, respec- 
tively. A scaling factor is introduced for pp collisions at 
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FIG. 5: Comparing pseudorapidity distributions among pp 
collisions at ^fs= 5500 and 200 GeV and 0-5% most central 
Pb+Pb and Au+Au collisions at ^JWn =5500 and 200 GeV. 
The panel (a), (b), (c), (d), (e), (f), and (g) is, respectively, 
for d, d, u, u, s, s, and g. In any of these panels the solid 
and open circles are for the Pb+Pb and Au+Au collisions, 
respectively, and the solid and dashed curves for the pp colli- 
sions. A scaling factor is introduced for pp collisions, see text 
for the detail. 



> io t 

u io [ 

O 10" 

sr 10 " 

E-h 10 



r% d 

i- (a) ^r^-~-^_ 

r| I I I | I I I | i'iTI I I +HH 


'(bj a"r 

iiiiiiMi>Tii7rrTTrir 


a 

r (c) v 


(d) u r 

iiMiiMiYri|TrrTTTir 


-1 " 1 1 H 1 | II T-| II +HH 

r s 

r 

PSa*.- ■ . . 


(f) s r 

r 

-^55*4 • • . J 


f <g) 


2 4 6 8 

Pt 



4 6 

Pt 



is because in the PACIAE model the nucleus-nucleus col- 
lision is decomposed into NN collisions and the NN col- 
lision decomposes further into parton-parton collisions. 
Meanwhile, we consider only 2 — > 2 parton-parton scat- 
tering processes and among these processes the elastic 
scattering is dominant. 

Similarly, one compares the different parton and anti- 
parton px distributions in pp collisions at y/s=5500 and 
200 GeV with the ones in Pb+Pb and Au+Au colli- 
sions at y / ij v ]v=5500 and 200 GeV in Fig. [6l respec- 
tively. However, in Fig. [6] we do not introduce any scal- 
ing factor. The same conclusion about the similarity 
discussed in the last paragraph may emerge if one exam- 
ines the shapes of parton and anti-parton pt distribu- 
tions in Fig. [HI Meanwhile, one sees in Fig. [5] that the 
discrepancy in slope parameter (apparent temperature) 
between pp and nucleus- nucleus collisions at the same 
cms energy decreases with cms energy increasing. 

In summary, we have used a parton and hadron cas- 
cade model, PACIAE, to investigate systematically the 
yield, (pseudo)rapidity distribution, and transverse mo- 
mentum distribution of partons and anti-partons in par- 
tonic matter assumed forming in the early stage of pp 
and nucleus-nucleus collisions at RHIC and LHC en- 
ergies. The discrimination, in the yields, 77 distribu- 
tions, and pt distributions, between different parton and 
anti-parton, RHIC and LHC energies, as well as pp and 
nucleus-nucleus collisions is discussed. We have turned 
out that the different parton and anti-parton approach 
to have similar spectrum with reaction energy increas- 
ing from RHIC to LHC. This is consistent with the fact 
of the (dressed) mass effect decreases with momentum 
(reaction energy/temperature) increasing. This result 
seems not supporting the assumption of different quark 
and anti-quark have the same pt distribution at RHIC 
energy. We have also argued that if the partonic matter 
forming in the early stage of Au+Au collisions at RHIC 
energy is sQGP, the one forming in the Pb+Pb collisions 
at LHC energy might approach fQGP. 

The transport model results of this paper seem to sup- 
port that the fQGP may properly have following physical 
features: 

• Gluons carry more than half fraction in total mo- 
mentum (energy). 



FIG. 6: The same as Fig. [5] but for transverse momentum 
distributions and no scaling factor is introduced. 

^=5500 and 200 GeV. That is 100 and 100 in panel 
(a), 30 and 150 in (b), 50 and 50 in (c), 30 and 150 in 
(d), 30 and 150 in (e), 30 and 150 in (f), and 30 and 
100 in (g). We see in Fig. [5] that, in globally speak- 
ing the shapes of parton and anti-parton 77 distributions 
in pp collisions are similar to the corresponding ones in 
nucleus-nucleus collisions at the same cms energy. That 



• Ratio of anti-parton to parton approaches one. 

• Flavor balance, such as 
u:d:s: c„. « 1 : 1 : j s : j c 
and uu : dd : ss, ... ~ 1 : 1 : 7 S 

(7s > 0.3 and j c > 10~ n for instance) is expected. 

These physical features might shed light on the searching 
for fQGP. 
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